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We present a complete solution of the soft x-ray edge problem within a field-theoretic 
approach based on the Wiener-Hopf infinite-time technique. We derive for the first time 
within this approach critical asymptotics of all the relevant quantities for the x-ray prob- 
lem as well as their nonuniversal prefactors. Thereby we obtain the most complete field- 
theoretic solution of the problem with a number of new experimentally relevant results. 
We make thorough comparison of the proposed Wiener-Hopf technique with other ap- 
proaches based on finite-time methods. It is proven that the Fredholm, finite-time solution 
converges smoothly to the Wiener-Hopf one and that the latter is stable with respect to 
perturbations in the long-time limit. Further on we disclose a wide interval of inter- 
mediate times showing quasicritical behavior deviating from the Wiener-Hopf one. The 
quasicritical behavior of the core-hole Green function is derived exactly from the Wiener- 
Hopf solution and the quasicritical exponent is shown to match the result of Nozieres and 
De Dominicis. The reasons for the quasicritical behavior and the way of a crossover to 
the infinite-time solution are expounded and the physical relevance of the Nozieres and 
' De Dominicis as well as of the Winer-Hopf results are discussed. 
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I. INTRODUCTION 

a 

S ' The problem of x-ray absorption and emission in metals has been widely studied experimentally as well 

1-^ . as theoretically in recent years. Various approaches have confirmed that both absorption and emission 

■ spectra display a sharp threshold reflecting the edge of the Fermi distribution of conduction electrons. 

O , Foundation £ipr a micrjoscopic description of the singula* edge behavior was laid by the pioneering papers 

O ■ of AndersonEJ, Maharn, and Nozieres and co-workersaQ. They demonstrated that the electron-electron 

Coulomb repulsion and a collective many-electron response to a sudden generation of a transient core 
hole - conduction electron pair by absorption (emission) of light are responsible for the existence of the 
x-ray edge singularity. Fundamentals of the threshold behavior are captured by a microscopic, tight- 
binding model proposed by Mahano and solved nonperturbatively for the first time by Nozieres and De 
DominicisQ It is now called the Mahan, Nozieres and De Dominicis (MND) model. 

There are two separable processes contributing to the singular edge behavior within the MND model. 
First, it is a rearrangement of the ground state of N noninteracting electrons jdue to a sudden change 
of a local, core-hole potential known as Anderson orthogonality catastropheB. Second, it is motion 
of an electron excited by absorption of light from a localized coreJevel to an extended state away from 
equilibrium and subjected to the so-called "final-state interactions" □. The calculation of soft x-ray spectra 
at the edge represents hence one of the simplest, real many-body problems allowing for an exact solution. 

Although the model MND Hamiltonian is rather simple and its exact equilibrium thermodynamics can 
be found in closed form for all temperatures and interaction strengths, the x-ray absorption and emission 
represent nonequilibrium processes to which sophisticated techniques have to be applied. Several methods 
have been proposed and developed in an effort to determine the edge structure of soft x-ray spectraQ. 
They can be divided into two classes: field-theoretic approach using many-body Green functions and 
Kubo formula on the one side and quantum-mechanical calculations based on the Fermi golden rule and 
transition amplitudes between one-electron states (Slater determinants) on the other side. 
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A breaJithrough in the description of the x-ray problem was achieved in the paper of Nozieres and De 
DominicisQ. They succeeded to separate the dynamics of the core and the band particles and thereby to 
reformulate the problem as an effective one-body theory. The calculation of a two-particle susceptibility 
needed for the absorption (emission) amplitude was reduced to evaluation of two one-particle Green 
functions. These two Green functions are related to the two separable processes contributing to the 
singular edge behavior of x-ray spectra. The rearrangement of the ground state of N noninteracting 
electrons due to a sudden change of a local, core-hole potential is described by the core-hole Green 
function. Motion of an electron excited by absorption of light is determined by a nonequilibrium Green 
function of conduction electrons. Separation of the dynamics of the local hole and the transient conduction 
electron is the bedrock of the approaches aiming at an exact solution of the critical threshold asymptotics 
of x-ray spectra. 

Quantum mechanical approaches treat only the conduction electrons as live particles. The local core 
states act as an external potential on which noninteracting conduction electrons scatter. A transition 
between the ground state of free electrons and ei|genstates of an electron gas iijua local potential determines 
the x-ray absorption and emission amplitudeaj. Combescot and Noziereaj proposed a determinagtal 
formulation in which the transition amplitude is a sum of determinants. Mahan and co-workeraZrt2l 
used analytic properties of the scattering T- matrix and represented the transition amplitude as a sum of 
dispersion integrals. The most comppeiheHsive approach using quantum mechanical transition amplitudes 
was achieved by Ohtaka and Tanabetll~E3la who used the Fermi golden rule and the Fredholm method to 
sum up a series of Slater determinants. These quantum mechanical approaches have a great advantage 
that they are applicable to the evaluation of the Anderson orthogonality catastrophe as well as to the 
final state interaction. The effect of bound states due to strong electron hole attraction can be estimated 
within this approach too. However, the final result is a sum of an infinite series the convergence of which 
cannot be strictly proven, at least with rigorous mathematical methods. Moreover the eigenstates of the 
electron gas in a local potential are not known explicitly and an approximate scheme to sum over the 
electron momenta using the formula of FriedelEj relating the eigenenergies with the phase shifts is usedEl. 
Since the quantum mechanical approach of Ohtaka and Tanabe is able to determine not only the critical 
exponents of the edge behavior but also the nonuniversal prefactors, it is now considered as the most 
complete and exact solution of the MND problem. 

Simultaneously to the quantum mechanical approach to the x-ray problem a field-theoretic scheme 
has been developed. Contrary to the former method the latter treats both types of particles of the 
MND problem as dynamical objects subjected to a dynamics governed by the MND Hamiltonian. The 
absorption and emission amplitudes are then represented via a Kubo formula using many-body Green 
functions. In analogy to other many-body problems such as Kondo, Anderson or Hubbard models, 
various field theoretic techniques has been applied. Weak-coupling perturbation theory was used in the 
early papers on the MND problemBfl. Firs.t nonperturbative solution within the field-theoretic approach 
was found by Nozieres and De DominicisQ. An asymptotic approximation of the authors and an exact 
solution of the resulting singular integral equation led to the confidence that the exact solution to the 
MND problem had been found. However, mathematical justification of the exactness of the applied 
asymptotic transformation is missing and only universal critical exponents can be obtained from it. The 
full exact solution would demand to calculate also the nonuniversal prefactors shaping the edge of the 
spectra. The construction of Nozieres and De Dominicis is up to now the only field theoretic solution 
capable to evaluate the critical exponents for both the Green functions in the x-ray problem. Later 
approaches concentrated mostly on the Anderson orthogonality catastrophe, i. e. on the evaluation of 
the Green function of the local quasi-mxticle. 

Rivier and SimanekEj and HamannEj circumvented the asymptotic transformation of Nozieres and De 
Dominicis and used the formalism of the Hilbert problem. They found nonperturbative critical exponents 
for the overlap between the Fermi sea and a new ground state of the electron gas in a local Coulomb 
potential. It became clear from these papers how a delicate problem it is to find the proper limit of 
infinite time (lattice volume). The long time asymptotics of the Green function in thCpMND and related 
Falicov-Kimball model was analyzed recently by means of the Wiener-Hopf method irO. There is up to 
now no analogous analysis of the nonequilibrium Green function of the conduction electrons reflecting the 
final-state interaction in the MND problem. The field-theoretic approach is considered at present as less 
complete than the quantum mechanical one. Moreover the critical exponent of the Green function of the 
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local, core electron (hole) calculated with the nonperturbative Hilbert-problcm (Wiener-Hopf ) trick does 
not match the quantum mechanical result. While quantum- mechanical and Nozieres and De Dominicis 
solutions predict a value (5(0)^/7r^, the Wiener-Hopf and Hilbert boundary-problem techniques produce a 
value (5(0)^/277^. Recently thje-author analyzed the difference between the Wiener-Hopf and the Nozieres 
and De Dominicis exponentala'E^ and concluded that the value of the critical exponent depends on the 
way how the effective lifetime of the transient electron-hole pair and the relaxation time are limited to 
infinity. 

In this paper we complete the field-theoretic solution of the x-ray problem based on the Wiener- 
Hopf trick and derive asymptotic limits for the Green function of the transient conduction electron and 
of the x-ray absorption amplitude. Thereby not only the critical exponents are obtained but also their 
nonuniversal prefactors revealing a number of new physically relevant features. Further on we discuss and 
explain in detail the differences between the finite-time and infinite-time approaches to description of the 
critical edge behavior. We show that the finite-time solution smoothly converges to the Wiener-Hopf one, 
but that there is a wide interval of intermediate times where deviations to the Wiener-Hopf asymptotics 
are apparent. We derive an exact formula for a long-time asymptotics for the Green function of the core 
hole including critical as well as quasicritical time regimes. We show that the quasicritical edge exponent 
is that of Nozieres and De Dominicis, i. e. twice the value of the Wiener-Hopf solution. The reasons for a 
quasicritical behavior and the way of a crossover to the infinite-time Wiener-Hopf solution are expounded 
and the physical relevance of the Nozieres and De Dominicis as well as of the Winer- Hopf approaches are 
discussed. 

The layout of the paper is as follows. Section H summarizes the basic mathematics of the x-ray 
problem, i. e. the fundamental integral equations and available methods for their solutions. In Sec. Ill 
we hint to an inconsistency between the critical exponents obtained from the orthogonality catastrophe 
and from the MND problem. In Sec. IV we present the full Wiener-Hopf long-time asymptotics for the 
relevant quantities in the x-ray problem. Section V analyzes the relevant long-time scales and shows 
the way how the critical or quasicritical behaviors with different exponents may arise and how various 
quantities reach their infinite-time limits. In section VI we discuss physical and practical consequences 
of the obtained results. 



II. MATHEMATICAL BACKGROUND OF THE X-RAY PROBLEM 

The aim of the paper is to discuss the source of the differences in the predictions of the critical 
exponents of the edge singularity. To avoid any ambiguities due to different definitions or formulas 
and to be able to discuss exact solutions we first summarize the basic mathematics of the problem 
to be studied. To determine the genuine long-time asymptotics of the x-ray absorption or emission 
amplitude, it is necessary to specify first the model within which we will work. We hence define the 
model and fundamental equations to be solved. Then we characterize mathematical methods for solving 
the fundamental equations that may lead to exact solutions, at least in particular limiting cases. 

A. Model and fundamental quantities for the x-ray problem 

The canonical model of the x-ray problem is that of Mahan, Nozieres and De Dominicis, the dynamics 
of which is governed by a tight-binding Hamiltonian 

Hmnd = ^k4ak - Ecb^'b -J^f^^Yl (1) 

k k,k' 

where b,b\ and ak,a]{j. are annihilation and creation operators of the core and conduction electrons, 
respectively. Although it is not necessary we confine the analysis only to the case of a contact potential 
Vkk' = y (lowest x-ray channel). Hamiltonian (|^) is a many-body formulation used in field-theoretic 
approaches. Since the local electrons have no kinetic energy they can be excluded from the problem. We 
then distinguish two Hamiltonians with conduction electrons only, the initial one 
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Hi = ^ ekaj^Ok (2a) 

k 

and the final one 

Hf ^^ewalak -Af-^^Vkk'alak^ . (2b) 

k k,k' 

The advantage of Hamiltonian (|l|) is that it equaUy describes the dynamics before and after absorption 
(emission) of light. 

Absorption of light is described by a transition amplitude obtained from the Fermi golden rule 

/•oo 

I{uj) = 2ReY, dte-*(-^^"^S~"+)*|(*/(7V + l)|T?|$(iV))|2, (3) 
/ ^0 

where $(iV) is the ground state of N noninteracting electrons from the conduction band, = Ep — Ec 
is the ground-state energy of the noninteracting system of N extended electrons and one local hole, and 
^' f{N + 1) are eigenstates of the Hamiltonian ( pb| ) with + 1 electrons labeled by the subscript /. The 

operator W = l/-^Sk ''^('^)'^k represents the process of creation of a new conduction electron due to 
absorption of light. Determinantal formulation uses this expression to evaluate the threshold behavior of 
x-ray spectra. 

In this paper we will work within the field-theoretic approach. Instead of (^ we use an equivalent 
Kubo formula with the conduction as well as with the local electrons 

^(^) = T72 E^^W^'lk')* / dt [e''^+*Xk',k(t) +e---*Xk.k'W] , (4) 
■'^ k,k' 

where Xk',k(i) = ^$(iV)|6T SMNDbHt)ak' {t)al{0)b{0) b^^{N)'j and uj± = uj ± irj. We denoted the 
S-matrix of the MND model in the interaction picture as Smnd — Tcxp|— i Jj^(ity(t)|, with the 

unperturbed Hamiltonian Hq — Hi — Ecb'^b. 

Fermi golden rule (^ uses quantum mechanical matrix elements that enable a direct access to the 
absorption amplitude. However the eigenstates of the Hamiltonian Hf are not known in closed form. On 
the other hand the Kubo formula (^) does not require knowledge of the eigenstates of Hf, but it demands 
Green functions to be known off the mass shell. The quantum-mechanical approach was exhaustively 
described inEl and we do not deal with it any longer. We will follow the quantum-field approach with the 
Kubo formula (^). _ 

Nozieres and De Dominici^ made a fundamental observation that the two-particle susceptibility 
factorizes due to the locality of the core particle into a product of two one-particle Green functions 



SMNDbHmo)] &t|cO(7V)\ U(N)\bT \Smnd iO,t) a^,{t)al{0)\ b^^N)) , (5) 



Xk%k(i) = 
<^>{N)\bT 

where Smnd (0, t) denotes the S'-matrix on a time interval [0, t\. The former function expresses an overlap 
between the ground states of Hi and Hj and is responsible for the orthogonality catastrophe. The latter 
one is a nonequilibrium Green function of the conduction electrons containing the final-state interactions. 
While the Green function of the local core hole was calculated by several field-theoretic techniques there 
is no nonperturbative, field-theoretic solution for the nonequilibrium Green function of the conduction 
electrons. We construct it in the next sections. 

We now determine fundamental dynamical equations for auxiliary Green functions from which both 
the averages on the right-hand side of (|^) can be determined. The average of the core particle can be 
represented as 
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igh{T) = (^^{N)\hT \SMNDb\T)h{Q)\ 6t|$(7V)^ ^exp{-z£;,r + C(r)} 



(6a) 



where C{T) is a sum of connected closed loops of the conduction electrons. It can be rewritten as an 
integral over the interaction strength and an auxiliary " open-line" Green function 



(6b) 



(6c) 



C{T)^ I dt' ( dWrit' ,t'^;\) . 
Jo Jo 

The Dyson equation for the Green function Txiti, t2] A) looks like@Ji3: 

TT{tiM-A) = Gc{ti^t2)-\ f dt'Gc{ti~t')TT{t'M-A) , 

Jo 

where Gc{t) = A^"^ Ek ^c(i, k) = (27r7V)~iEk/ duje^"^* [lu - + ir]sgnuj]^^ is the diagonal (local) 
element of the causal Green function of the conduction electrons, A is an intermediate interaction (Ae 

The Dyson equation for the nonequilibrium Green function of the conduction electrons has a structure 
analogous to (gg). We use a representation 



(7a) 



(7b) 



rT(T,0;k',k)= -^($(7V)|fer \Smnd (0, T) ak'(?^)4(0)J b^MN)] 

= Gc(T, k)(5(k - k') - V" / dt'Gc{T - t', k')TTit', 0; k) , 

Jo 

where the function rT(ii,i2;k), T > ^1,^2 > fulfills a closed equation 

rT(ii, ^2; k) = Gcih - t2,-k)-V [ dt'G.ih - t')TTit\ ^2; k) . 

Jo 

Equations ( |6c| ) and (^) are the desired exact Dyson equations, to be solved. They are identical with 
equations (17) and (22) of the Nozieres and De Dominicis papera. The kernel of these integral equations 
is the diagonal element of the Green function of noninteracting conduction electrons and is explicitly 
known contrary to other many-body problems. The only evidence of many-body effects here are the 
closed loops in the core-hole Green function and the transient (nonequilibrium) form of the band-electron 
propagator. 



B. Basic techniques for solving Dyson equations in the MND problem 

The above field-theoretic approach reduced the MND problem to two separate integral equations with 
the same integral kernel. Exact solution of the x-ray problem then equals an exact solution of equations 
( |6c|) and (|7b| ). They can be written in a generic form 

TT{h,t2;n)=T^°\h-t2;n)-X [ dt'G^h - t')TT{t' ,t2; n) , (8) 

Jo 

where A £ [Opk"lis interaction strength and k is a parameter depending on the type of the Green function 
to deteHH-inellMcy. For T < cx) it is a Fredholm integral equation and can be solved using the Fredholm 
methodB. With it we can find a solution in the Hilbert space of square integrable functions on the interval 
[—T,T], L2{—T,T). To obtain a unique solution to (||) we have to fulfill the following two conditions: 

Fl: Integral kernel Gc(ii — ^2) is a Hilbert- Schmidt operator in variables ti,t2, i. e. 

\\Gc\\hs-= [ dtidt2\Gc{ti~t2)\^ =T/2 [ dt\Gc{t)\^ < 00 (9) 

Jo J-T 

and 
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F2: Inverse interaction strength —A ^ is not a characteristic number of the integral kerneL i. e. it is 
not from the discrete spectrum of the operator Gc defined through its matrix elements 

{ti\Gc\t2) Gc{ti — 12). It means that = —X^^ip has no nontrivial solution. 



Then the solution to (0) can be written as a spectral representatioi 



J2^ 



rF{ti,t2;K) ^ 1 (10) 

i + AA„ 

where $„ are eigenfunctions of G with eigenvalues A^^. The brackets ( ) stand for the scalar product 
on L2{~T,T). The eigenvectors \^n) form a complete set, i. e. J2n\^ri){^7i\ = 1 on L2{-T,T) and 
= ^n^n^\^n){^n\ in average On L2{—T,T). Both eigenvalues A~^ and eigenfunctions depend on 
the interval length T. For finite time intervals we can use the Fredholm determinants and their analytic 
properties to continue the results of perturbation (Neumann) series to arbitrary interaction strength. 
However, if the interval length T tends to infinity the spectrum of the integral kernel Gc(ii — t2) goes 
over into continuum. Hence to keep the interaction strength out of the spectrum it is necessary to add a 
small imaginary part irj to it, i. e. A ^ A + irj. Otherwise it is impossible to perform the limit T — )■ 00. 
Mathematically this imaginary part secures the convergence of the formal limit T 00. To keep the 
theory unitary it is necessary to push 77 — > in the end of the calculations. Physically a small imaginary 
part in the interaction reflects an effective half-time of the exponential decay of the Coulomb potential 
between the core hole and conduction electrons defined in our units, {h = 1), as r = l/rj. The time r — >■ 00 
defines a new "infinite" time scale, effective lifetime of the transient core hole (conduction electron). It is 
an experimental quantity depending on the recombination probability in the x-ray process. It is essential 
for the Fredholm method that although the relaxation time T 00, the ratio T/r — > 0. Most of the 
existing solutions to the MND problem are based or equivalent to the Fredholm method characterized by 
the ratio T/r 0. 

Another technique how to solve the MND problem, i. e. the integral equation (H) in the asymptotic 
limit r — > cx) is the Wiener-Hopf trick. If the interaction strength — A becomes a characteristic number of 
(H) and the spectrum of the integral kernel Gc{t) is continuous, then (^) and (|l^) become divergent. We 
have to replace the Fredholm with the Wiener-Hopf method based on the Hilbert boundary problemE3. 
The Wiener-Hopf technique solves (^) exactly and at weak coupling uniquely under two assumptions: 

WHl: i2(~T, T) = L^ ® L+ := L2(-r, 0) © ^2(0, T); 

WH2: if /,g e L+{L-) then (/ * g){t) := J^^dt'f(t - t')g{t') is orthogonal to L-{L+). 

These conditions are strictly fulfilled only in the limit T = 00. At T < 00 the spectrum of Gc{t) is 
discrete and WH2 is violated. The spectral representation (|l^) is in the Wiener-Hopf solution replaced 
with another decomposition 

~ fW(u;2;^) ^-{^2) , 

r+{uji,uj2;K) ^ — , iia) 

i(cj2 - uji ~ irj) <I'+(wi) 
where the "spectral sum" (integral) appears in the exponential functions 

^±{x) - exp -— / dy — ln[l + AG,(j/)] . (lib) 

Here G^{iu) := J^^dte^^'e^'^^'^Gcit) and 

r±(wi,W2;K) := /^diid<2^(±ti)e*"i*ie-*"2*2e"''(l*il+l*2l)r(ii,i2;'«) are Fourier transforms. Contrary 
to the Fredholm method the Wiener-Hopf trick does not anticipate knowledge of the eigenstates and 
eigenvalues of the integral kernel. The small imaginary part rj has the same meaning as in the Fredholm, 
finite-time approach. It secures the convergence of the Fourier integrals used in the Wiener-Hopf solution. 
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However, the parameter 77 has quite a different relation to the relaxation time. The relaxation time T is 
strictly infinity in formulas ( |ll|) whereby the parameter 77 must remain finite. We hence have rjT — > 00. 

Leading long-time asymptotics T 00 within the Wiener-Hopf solution is obtained by replacing 
the infinite time interval by a finite one in the solution (pi]). I. e. the Fourier integrals come over into 
Fourier series-aad we obtain in leading order of for the Matsubara frequencies w„ = {2n + l)TrT and 

r+( -A-^A r(")(^„;^) 

4 1 + XGc[Vm) 

r+l \ -1 f(")(^„;k) $+(c.„) 

r;^ z/„,w„;k) = — — -— -, 12b) 

where is a projection of the Wiener-Hopf solution F+ onto the interval [0, T]. The exponential function 
$+ on a finite interval splits into two, 

-In 1 + AG,(;.„0) -7^ In ( 1 + AGc(c^„) ) \ (13) 

2 V ) T ^ i(y^- i^n) V /J 

and <i>+(aj„), where fermionic and bosonic frequencies interchange their roles. 



III. INCONSISTENCY IN THE EXISTING INTERPRETATION OF THE RESULTS FOR 

THE CRITICAL EDGE EXPONENTS 

It is generally accepted that the various existing approaches to the MND problem confirm the pio- 
neering result of Nozieres and De Dominicis on the critical edge exponentsu. We now show that it is 
not quite t he . ca se and that there is a discrepancy if the results obtained for the AjQdetsmi orthogonality 
catastropheli3ii3 and for the critical behavior of the Green function of the core holeO'Enid'Ej are compared. 
The discrepancy between the two existing results becomes apparent if we systematically use only one and 
the same construction scheme on both quantities. 

In the orthogonality catastrophe one has to evaluate a matrix element (<I>(iV)l\['Ar(t)) in the limit 

t ^ 00. Here |^jv(^)} = cxp ^—HImnd^ \^ i^))- ^ simple calculation yields 



$ (N) I'fNit)) = e-*^-* (N) \SMND{t, 0)1$ (iV); 
= e-*^^*^$(7V)|rexp|-zjJ^*di'T>/(<')| l$(^)^ = -ie^'^^'Ghit) . (14) 



The last equality follows from (6a) and the fact that the core electron does not propagate in the MND 
model. From the other side, if we start from the Fermi golden rule (^ and consider only the vacuum 
contribution, we obtain 



^e-*(^^-^«)*|(*/(iV)l$(Ar))|2 =^e-*^^*($(iV)|exp{iiJoi} |*/(iV)) 
/ / 

x{^fiN)\e^p[-iHMNDt} mN)) = e-'^=*($(iV)| exp {iHot} exp [-iHMNot} MN)) 



= e-'^^*mN)\SMND{t,0)\<i>{N)) - -iGnit). (15) 

It hence means that the overlap calculated in the orthogonality catastrophe ( |l^ is identical with the 
Green function of the local, core hole used in the Nozieres and De Dominicis and other field-theoretic 
approaches starting with Kubo formula (jj). Now the calculatioa. of the transition amplitude in the 
orthogonality catastrophe based on the Hilbert boundary problemlla leads to the following asymptotics 
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|($ (iV)|*Ar(t))| ->Cexp 



27r2 



(16a) 



with (5(0) = Imln 1 + V^G(— i??) > 0. The result of Nozieres and De Dominicis and othera3 for the 
Green function of the core hole reads 



\gH{t)\^Cexp 



(16b) 



If we accept the resuh (KW) as correct, we then get in conflict with the critical exponent of the Anderson 
orthogonality catastrophdH. If we accept the result for the overlap amplitude we do not recover the result 
(16b), since ( p^ ) and (15|) are simple exact relations derived without any special assumptions. 

lLis.|Worth commenting on the determinantal approach using the eigenstates of the final Hamiltonian 
(2b)Llol3. Their construction is based on two assumptions. First, they assume that ($ (N) |5'Ar(t)} — > 

{N)\'^'Q{N)) if i ^ oo, where |vI'o(iV)} is the ground state of the final Hamiltonian Hf with N 
conduction electrons. It means that adiabatic theorem is assumed to be valid. Second, it is shown that 



each term in the sum over the eigenstates of Hf from (|15|) contains a factor | ($ (N) |^o(-^)} | . It is tacitly 
assumed (but not proven) that the infinite sum converges in the infinite-time limit and hence the square 
of the overlap between the two ground states appears also in the Green function of the core hole. This, 
however, would be in conflict with the equivalence, of the Fermi golden rule (^ and the Kubo formula 
(^) proven in various books on condensed mattered. The explanation why the determinantal approach 
failed to disclose the discrepancy between the orthogonality catastrophe and the x-ray results is that the 
second assumption about the convergence of the sum over the eigenstates of Hf in the long-time limit 
does not hold. The interchange of the summation and the limit to infinite times is not justified. We see 
that our present knowledge of peculiarities of the edge behavior is insufficient and the problem deserves 
further analysis. 



IV. LONG-TIME ASYMPTOTICS USING THE WIENER-HOPF TECHNIQUE 

In this section we show how to use the Wiener-Hopf solution from the infinite-time limit ( p|) to 
evaluate long-time asymptotics of the relevant quantities in the MND problem. They are the Green 
function of the core hole, the nonequilibrium Green function of the conduction electrons, and the x-ray 
amplitude. Since the Green functions arc singular only at zero temperature, we give an explicit derivation 
for this case. 

As we already mentioned, the long-time asymptotics is derived within the Wiener-Hopf method if 
we take the infinite-time solution ( |ll| ) and project it onto a finite interval [— T, T]. The functions are 
periodically continued beyond this interval. Strictly speaking the Wiener-Hopf formulas (p^), ( p^ are no 
longer exact on a finite interval, but we show in the next section that the exact solution deviates from the 
Wiener-Hopf one only in higher orders of the inverse length of the interval, T^^ and that the finite-time 
solution converges to the Wiener-Hopf one at sufficiently long times. 



A. Green function of the core hole 



Green function of the core hole was extensively studied in the framework of the Falicov-Kimball model 
in Ref.EII. The Falicov-Kimball model differs from the MND one only in that the impurities are densely 
distributed on the lattice. The method of solution is otherwise the same. Representation (|6a| ) is used 
together with the auxiliary Green function Txit' , t'_^_; A) and defining equations ( |6b[ ) and (|6c[). These are 
exact equations that, however, do not possess exact analytic solution in the whole (infinite) time range. 
The long-time asymptotics calculated with the Wiener-Hopf solution has the following form 
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--y- 



$+K) 1 + XG,{V^) $+(i^m) 1 + AGcK) 



(17) 



Since we are interested only in leading-order long-time behavior, we replace the discrete sums with 
continuous integrals whenever possible, i. e. when thajntegrals converge. The discrete sum is a regular- 
ization needed to make the infinite-time limit scnsiblec£l. First of all the exponential function <i>+ can be 
represented as 



$+((jj) — cxp < — 



duj' 



-Iniln 



TT J UJ ~ uj — irj 




1 + \G{uj' + if]) 



and the sums in ( |l7|) can be rewritten into 



C{T)^— [ dLo5{uo) 



dX—P 
2-K 



(18) 



do; 

1 . 1^ 



(19) 



where I+{ui) := 0{tj){I{uj) - /(0+)) + 9{-uj) {I{uj) - /(O^)) and S{uj) = Imln 1 + VG{uj - irj) 

G{uj' + if)) 



with 



<iw'$+(w') 



Im- 



1 + XG{uj' + iri) 



(20) 



The first term on the left-hand side of ( [l9| ) is a shift of the gijound-state energy of the conduction electrons 
due to the core-hole potential matching the Fumi theoremn3. The second one is a logarithmically singular 
contribution causing a critical algebraic long-time decay of the Green function of the core hole. The 
logarithmic divergence of the integral with the function /+(w) appears due to a jump in the derivative of 
this function at the origin. When evaluating this jump we obtainllZI 



u 

^(^) = V / duj5{u)-^5{Qf\niT, 



(21) 



where ^ is an effective bandwidth of the conduction electrons. p^his result is in accord with the calculation 
of Hamann done for the Anderson orthogonality catastrophet3. 



B. Nonequilibrium Green function of the conduction electrons 

Although there are various approaches to the orthogonality catastrophe or the Green function of 
the core hole, there are only a few techniques to calculate the nonequilibrium Green function of the 
conduction electrons (0). It is surprising but a full, field-theoretic solutiorufor this function has not yet 
been published. Except for the series of papers by Nozieres and co-workersQ'Q the solution for this Green 
function are based on the quantum-mechanical golden rule or Slater determinants. We present now a 
field theoretic, Wiener-Hopf solution to this nonequilibrium function in the long-time limit. 

To reach the long-time asymptotics with the Wiener-Hopf method we must solve equation ( [7b| ) for 
finite time intervals, then put ti —T and t2 = and only after this perform the limiting process T ^ oo. 
The Wiener-Hopf solution as given in (111) applied to eq. (7b) reads 
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r+(i/m,a;„;k) 



1 + VGcil^m) 

-1 Ge(a;n;k) $+(a;„) 



(22a) 
(22b) 



with <i>+(^'„i) and $-|_(a;„) defined in (|l3|). These expressions are regular if we properly define the meaning 
of the Green functions G'c(j^m = 0;k) and Gdv = 0). These functions are not uniquely defined due to 
poles (cut) of these functions at the real axis. To avoid this difficulty we redefine any function with zero 
Matsubara bosonic frequency as 



/(^,„ = 0)^-[/(0+) + /(0-)] 



(23) 



The Wiener-Hopf solution to the full nonequilibrium Green function of the conduction electrons ( |7a| ) 
reads 

TriT, 0; ki, k^) ^ e"^'"''GT(z/™; ki, ka) - ^ ^ e"^"''GT(c.„; ki, k2) 



l + VGc{LJn) *+K0 



+ {Vrn ^ Wn) 



(24) 



Here we denoted the equilibrium Green function of the conduction electrons with the core-hole potential 
V as 



Gt(c^; ki,k2) = -Gc(c^;ki) 



(5(ki - ka) - 



V 



1 + VGciuj) 



Gciio-M) 



(25) 



Relations ( p^ - (|2^) complete the Wiener-Hopf solution for the nonequilibrium Green function of the 
conduction electrons in the MND problem with a contact core-hole potential. It is worth noting that this 
solution is unique for a sufficiently weak interaction V, i. e. when no bound states interfere. 

It is easy to confirm that Foo^oo, 0; ki, k2), obtained from ( p^ in the limit T — > oo, vanishes, since 
the fermionic and bosonic contributions exactly cancel each other. Only the terms proportional to the 
powers and higher do not compensate and render a finite contribution. They contain singular terms 
dominant in the long time limit. To obtain the most singular expressions in the long time limit, it is 
sufficient to consider only contributions from the lowest, both fermionic (n = ±1) and bosonic (m = 0), 
frequencies with the function in the denominator. Namely only then we obtain a logarithmically 
divergent contribution to the sum. The contribution from m = (c.f. Eq. (|2^)) is one half of the sum 
over the lowest fermionic frequencies with n = ±1, while the higher bosonic frequencies m = ±1,2, ... are 
exactly compensated in the long-time limit by the fermionic ones n = ±2,3, respectively. 

We assess the long-time limit of the sum in the exponent of the function $_|- from (^2|). The sum 
turns to a principal value integral. The causal Green function Gdi^) has a jump in its imaginary part 
at = where the integrand linearly diverges. These two coincidences lead to a logarithmic singularity 
in the limit T — > oo. A straightforward calculation yields an asymptotic form of the functions 

$+(tj = ±vr/T) <^>+{i^ 



0^ 



exp | — In |^T||, where ^ is a suitable (nonuniversal) 



cutoff for the asymptotics lu We assumed 77T 00, the Wiener-Hopf regime. Inserting these 

singular contributions to (E4) we obtain 



rT(r,0;ki,k2) = — (^T) TT Re 



V 



P 



duo $_(cj) 
27ri LO 



Gd^-M) 



(26) 



as leading, singular term in the limit T ~* 00. We used 
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The critical exponent is again one half of that of Nozieres and De Dominicis and other finite-time (de- 
terminantal) approaches. 

Expression ( p6| ) can still be simplified. The integral with the fmiction $_ can explicitly be evaluated 
using analytical properties of and Gc- We have 



P I ^^Im 



1 



w - + iri 



-1>-(ek.). 



(28) 



Using this result we can write a final expression 



5(0) 

rT(T,0;ki,k2) = -— (CT) 7T 



V 



(ekj^-(£kj 

Ek, 



Re 



„i5(0)/2 



fki + «?7 



(29) 



Note that momenta ki and k2 do not enter the singular part of the nonequilibrium Green function 
symmetrically. The reason for this lies in asymmetric treatment of the ends of the integration domain 
in the fundamental integral equation (7b). While at t = the system undergoes a sudden change, i. e. 
the core potential is immediately switched on, the other end is pushed to infinity, T — > cx), and the 
core potential is being adiabatically damped by a small exponential factor e~''-^ to make the long-time 
limit meaningful. There is no abrupt change in the long-time tail. A sharp edge (represented by a step 
function) and a long-time tail dominated by the long-time asymptotics of the integral kernel Gc{T) are 
equally important for the existence of the edge singularity. If the sharp edge at i = were removed, 
the function rr(T, 0; ki, k2) would reduce to the equilibrium function GtCT; ki, k2). If the integration 
interval were kept finite we could not get to the Fermi energy closer than to a distance ttT and no 
logarithmic divergences would show up. 



C. Absorption amplitude 

Having derived expressions for both the Green functions we can now construct the absorption ampli- 
tude (|^). We first write down the two-particle susceptibility Xk,k'(^) i'^ the asymptotic limit of long times. 
Using formulas (^), (||), ([21}), and (^6|) we explicitly obtain for leading (singular) long-time asymptotics 

u\ ^"^^ rcrr^a f ..p 6l(ek')*-(ek')T, 
Xk,k'(i) = ^^pjT (CT) e^p{~i{Ec- £)T} Re 

Zl Ek' 

where £ = J_^duj5{uj) and a = 5{0)/tt - 5{'df/2TT'^. The factor i appears due to the definition of 
the Green function of the core hole (|6a|). This expression is exact (within the Wiener-Hopf infinite- 
time approach) not only for the critical exponent of the algebraic decay but also for the nonuniversal 
prefactor, the time- independent part of (^0|). To obtain the absorption amplitude I {to) we need to know 
the susceptibility x ^lot only in the long-time limit but on the whole positive frequency axis. However, if 
the critical exponent a from ( |30| ) is positive, i. e. the absorption amplitude diverges at the Fermi level, 
then the long-time behavior of the susceptibility also determines the edge behavior of the absorption 
amplitude. In the isotropic case with contact core potential we can derive an exact expression for the 
absorption amplitude in the vicinity of the edge. We insert (|3^) into (^) and realize that from the real 
part in (^0|) only the Fermi energy (e = 0) survives. Then the absorption amplitude is 



Ck + ir] 



(30) 
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I sing (l^) = 2Im 



2 ./o £ Jto t 



(31) 



where p(e) is the DOS of the conduction electrons and w(e) — 'w(e)* = w(ek). We stress that Iging 
is only the singular part of the absorption amplitude, i. e. a > 0. We used a new parameter Tq as 
a lower bound for the time integration, since only there the long-time asymptotics of the electron-hole 
susceptibility x dominates. The value of this parameter depends on the range of applicability of the 
Wiener-Hopf solution. We will discuss the problem in the next section. At present we know that at least 
To « l/ry. In the Wiener-Hopf regime then Tq/T and AujTq < 1, with Auj ^ ui - + £■ We can 
extend the time integration onto the whole positive axis and perform the integral explicitly. It yields 
exp{ia7r/2} [C/ {uj-Ec + E)^ T[a). 

The absorption amplitude must be positive, otherwise it is unphysical. It is the expected, 
only for ijj > Ec — E > Q ■ We explicitly have 



(5(0) an f°° $ (e) 

h^ng{oj) = 27rFp(0) sin(^) sin(il-)r(a)«;(0) / dew{e)p{e)—^ 







A, 



(32) 



The integral in (^2|) is finite, i. e. does not bring any new divergence into the absorption amplitude. 
The critical edge exponent in the Wiener-Hopf solution is one half of the critical exponent obtained from 
the Nozieres and De Dominicis and finite-time approaches. Formula (^2|) is the Wiener-Hopf expression 
for the asymptotic limit w — > of the absorption amplitude provided the leading low-frequency term is 
singular, i. e. a> 0, and the frequency is sufficiently close to the threshold, i. e. AwTq ^ 1. 



We can make a few general statements about the threshold behavior of the absorption amplitude (32). 
First of all the singular part of the absorption amplitude is very weak at weak coupling, namely it is 
proportional to {V p{Q))^ . Second, initial and final states are differently treated in the absorption process. 
The initial core electron can be excited to any empty state (e > 0) whereby the energy is not conserved 
due to instantaneous absorption. The excited state decays in the long-time limit towards a final state 
where the conservation of energy is restored and the additive conduction electron settles onto the Fermi 
surface (e = 0). A new ground state is reached. This is a very different result from the conjecture that the 
diagonal ij>aijt of the two-particle susceptibility, Xkk, yields the principal contribution to the absorption 
amplitudeB'El, i. e. the initial and finite states enter the absorption amplitude symmetrically. 



V. RELEVANT TIME SCALES IN THE X-RAY PROBLEM AND TRANSITION FROM 
FREDHOLM TO WIENER-HOPF SOLUTIONS 

In the preceding section we derived the Wiener-Hopf long-time solution for the relevant Green functions 
of the MND problem and also the edge behavior of the absorption amplitude. We obtained a critical 
behavior with a power law but with critical exponents that equal one half the values derived by Nozieres 
and De Dominicis and most other approaches. To understand the difference it is necessary to distinguish 
carefully relevant time scales in the x-ray problem and to determine in which limiting cases either solution 
is applicable. 

The MND Hamiltonian (|l|) contains two relevant finite energy or time scales. It is an effective band- 
width ^ related to the kinetic energy and the core hole - conduction electron interaction strength V . In 
our notation V > Q. Apart from this model-dependent energy scales there are a few "infinite-time" scales 
refiecting experimental realization of the process. First of all it is the relaxation time T determined by 
the waiting time at which a new equilibrium after the absorption of light is reached. This time determines 
the length of the time interval on which we solve the fundamental Dyson equation (^). There is also an 
eff'ective life-time of the transient pair of a core hole and a conduction electron, r = I/77. Apart from 
these two infinite times there is also an infinite scale connected with volume of the sample, or with the 
number of lattice sites N . The last very large scale is the inverse temperature (3. All the scales, not 
defined from the Hamiltonian, must be limited to infinity to reach an ideal critical behavior, the edge 
singularity. The actual critical asymptotics may be "unstable" and may depend on the order in which the 
limiting infinite values are reached. We compare only the two general approaches discussed in Sec. II. B, 
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i. e. the Fredholm and the Wiener-Hopf one, since they show the differences in the criticai exponents we 
want to understand. 

We k;now from the discussion in Sec. II. B that the Wiener-Hopf solution is exact if t/T ^ 0, i. e. if 
the relaxation time is much larger than the effective life-time of the core hole. This is, however, only a 
sufficient condition for the Wiener-Hopf infinite-time regime to set on. We would like to find necessary 
conditions distinguishing the Wiener-Hopf solution from the Fredholm one. 



A. Small parameters for the Fredholm and Wiener-Hopf solutions of the Dyson equation 

Restrictions on applicability of the Fredholm solution to the equation (||) in the long-time limit 
come from the assumptions Fl and F2 of the Fredholm method. The Fredholm method is based on 
a (Neumann) perturbation expansion that must have a finite convergence radius. Hence using we 
obtain an upper bound 

= V^\\G\\jis =T/2 [ dt\Gc{t)\^ ^n^V^p{Of^T < M , (33) 

J-T 

where M is a constant of order unity. Further on the spectrum of the kernel of a Fredholm integral 
equation must be discrete. It means that the distance between the nearest energy levels must be greater 
than the resolution of the measuring instruments. Combining this result with the restriction due to the 
effective life-time of the core hole we obtain 

T < JV^/^ , r</3, T«r=-. (34) 

V 

The above conditions are sufficient to secure the Fredholm solution is reliable and valid. However bound 
( p^ ) is very restrictive in long times. It tells us that we cannot reach the infinite-time limit for a fixed 
interaction strength V. We can assume A from (|3^) to be a small parameter for the Fredholm solution. 
It means that the smaller A the better the Fredholm solution reproduces an exact one. The latter three 
bounds say that the edge exponent for the absorption amplitude governing the infinite- volume limit 
A/" ^ oo must be that of the Wiener-Hopf solution. The same holds also for the limit ^ cxd for Au; « rj. 

It is more elaborate to find a small parameter for the Wiener-Hopf solution, i. e. an analogue to 
A from (|3^). The Wiener-Hopf trick leads to a nonperturbative solution and the small parameter will 
not depend on the interaction strength explicitly. Deviations of an exact solution on a finite interval 
appear due to violation of the assumption WH2 from Sec. II. B. A small parameter for the Wiener-Hopf 
solution is a measure of these deviations. The problem with the application of the Wiener-Hopf trick 
onto a finite- interval integral equation lies in different decompositions of Fourier transforms diagonalizing 
convolution. For a semi-infinite interval we decompose any function into its positive and negative parts, 
i. e. functions the Fourier transforms of which are analytic in the upper (lower) half-planes of complex 
energies, respectively. For a finite interval we have to decompose a function into its fermionic and 
bosonic parts in order to factorize (diagonalize) convolution in time. We relate the two decompositions 
and quantify deviations from the Wiener-Hopf solution in the long-time limit. 

Let us have a function F{t) G L2(— oo, oo). Let F± be its projections onto positive and negative time 
half-axis, respectively. We now cut off the function F at ±T and continue the function periodically onto 
the rest of the real axis. We denote such a function as F{T; t). It is sufficient to consider the argument 
t E (— r, T]. We use a variable C = ±1 for bosonic and fermionic projections, respectively. Then 

F+{T;t) = [F+{T;t) + Ce-^^e{t)F^{T;-T + t)] , (35a) 
F^- [T- 1) = [F_ (T; t) + Ce-'^^9i-t)F+ (T; T + t)] , (35b) 

where again t G (— T, T]. Note that the functions F^{T;t) have two contributions. Each is centered 
around the either edge of the interval [0, T]. The Wiener-Hopf solution does not treat the edges symmet- 
rically and the contribution from the upper edge is suppressed by a small exponential term e~^'^ . Using 
(|35| ) we define a general difference function 
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AF(r; t) = e"''^ [Oit)F^ {T;-T + t)+ e{-t)F+ (T; T + t)] 



(36) 



useful for building up a small parameter for an expansion around the Wiener-Hopf solution. 

The Fourier transform of the above functions and its analytic structure in the plane of complex energies 
is essential for the Wiener-Hopf solution. Using the fact that -^'±(2:) are analytic in the upper (lower) 
half-plane and supposing that they vanish at the respective infinity, we obtain from (|35|) 



F±(T;c.) = i[l + Ce-''^T*-^]F±(a;), 



(37) 



where F± are Fourier transforms on the infinite interval (—00, 00). The difference function, when specified 
to the integral kernel of the interal equation (0) has a Fourier transform 



'G„(w) + e~''^^G+(tj) 



(38) 



where the projections G± are defined analogously to 



and 



duo'- 



UJ' 



UJ ■ 



IT] 



-Ge(w') = - 



cLu 



,ImG(w' + irj) 



G_(c^) =p_[G,](c^) 



1 

2Td 



duj'— Gc(w') 

, LU — UJ + irj 



duj 



bj' — LO — ITj 

,Im G{uj' + irf) 



oj' — UJ + irj 



(39a) 
(39b) 



To find a small parameter for an expansion around the Wiener-Hopf solution we have to consider the 
fundamental Dyson equation (^-jn frequency representation. The equation splits into a couple mixing 
the bosonic and fermionic partsEZI. We can formally write the pair of equations as 



^(0) 



2VG,jT+ 



where, if we assume only the active variable. 



^(0) 



(40) 



(w„) ^pf[T]{uj„ 



m— — 00 ^ ^ 



(41) 



are projections onto positive (negative) time axis within the interval [— T, T\. From the solution (|4y) again 
only the positive projection has physical meaning. We see that the bosonic and fermionic functions 
are mixed. The difference between bosonic and fermionic functions in the plane of complex energies is 
given by the difference function (|3^). The bosonic function vanishes at fermionic frequencies a,nd vice 
versa. If the functions r&(ti-') and r/(w) were identical, the Wiener-Hopf solution would be exacttZl. 
In the long-time limit (in leading order) the projectors can again be expressed as integrals 



G 



duj'- 



UJ' 



bj ^ irj 



G, 



{uj) + Cp- 



AG7 



Gc{cj') - C^Griuj') + CAGt(w) 



(42) 



The difference between ( ^2| ) and the infinite-time expression is only in the function AGt{i^)- 

As we already mentioned the Wiener-Hopf solution on finite intervals breaks down because the as- 
sumption WH2 is violated. It means that convolution does not remain within positive or negative 
subspaces Lf- A measure of this "non-orthogonality" can be defined as 



G, 



G, 



(43) 



It is evident that k^(w) = at T = 00. The bosonic and fermionic functions (labeled with subscript ^) 
loose their meaning when going over to infinite time. We have to come to positive and negative projections 
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on the whole real axis. We sum over the index C to obtain the Wiener-Hopf solution r_|_ independent of (. 
We can introduce a new C- independent small parameter kt — '*C' ^fter a few simple manipulations 
we obtain 



1 



duj'- 



P- 



Lu' — LU + irj 



(44) 



We use explicit representation ( p3q ) for the function AGxii^) and perform the integral over x analytically. 
In the remaining expression we substitute x — loT and realize that only small values of x contribute in 
the long-time limit. The contributions from higher frequencies are canceled due to negative interference. 
If ^ is an effective bandwidth we finally end up with 



KT{io) = p{Of (In^T) 



2 e 



-2jjT 



dx 



smx 

X — ujT 



(45) 



where p(0) — — l/7rlm Giirj). Since we are interested in the limit T —^ oo with a fixed frequency to, we 
assume that luT ^ 3> 1. The integral is independent of frequency and can be bounded from above by 



46-^"^ 
7r2(eT)2 



(46) 



A dimensionless small parameter is obtained if we devide the norm of kt{^) on an interval [—CT,S,T] by 
the norm of the function G-\-{lu)G^{lu) on the same interval. The desired small parameter then is 



A 



(47) 



We see that the finite-time solution approaches the Wiener-Hopf one r_|_ algebraically and the small 
parameter is proportional to (^T)^^. The Winer-Hopf solution is hence stable in the long-time limit and 
we can systematically expand around it. 

Having found small parameters for both the Fredholm and Wiener-Hopf solutions we can distinguish 
their domains of validity. Comparing ( ^3| ) and ( pl| ) we see that the Fredholm, perturbative solution is 
applicable for ^ 1 while the Wiener-Hopf one in the opposite limit ^ 1 . There is hence a (smooth) 
transition from the Fredholm to the Wiener-Hopf solution if both the small parameters are approximately 
equal. Comparison of the parameters A and A as obtained in ( |33| ) and (|4^) yields an estimate for the 
transition time Tc and frequency lUc 



Atj, 




7r4T/2p(0)2 



(48a) 
(48b) 



The long-time Wiener-Hopf solution to the Dyson equation (||) becomes asymptotically exact if the 
relaxation time is sufficiently larger than the crossover value Tc or when the frequency is correspondingly 
smaller than Awr. 



B. Critical edge exponents 



We determined the domains of validity of the Fredholm (short-time) and the Wiener-Hopf (long-time) 
solutions to the equation (||). They are opposite edges of the parameters and rjT. The former is 
characterized by £,T <C 1 and rjT ^ 1, while the latter is valid if ^ 1 and r/T ^ 1. The two methods 
are hence, in some sense, complementary. Formally it is possible to continue weak-coupling perturbation 
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theory term by term to the long-time regime and one obtains an estimate for the critical edge behavior. 
This estimate coincides with the weak-coupling expansion of the solution of Nozieres and De Dominicis. 
However, perturbation theory in the long-time regime is very sensitive to the order in which the earlier 
discussed various large scales are limited to infinity. One has to be careful in using perturbation expansion 
in the long-time limit. On the other hand the stability of the Wiener-Hopf solution does not yet mean that 
all the long-time asymptotic results are determined from the Wiener-Hopf solution r_|_ only. To derive 
the algebraic decay of the Fredholm solution to the Wiener-Hopf one, we summed over the subscript 
Such a symmetric sum is well defined at infinite time even without a damping factor. The antisymmetric 
part, namely the difference between the bosonic and fermionic functions, rapidly fluctuates between r_|_ 
and and allows for an infinite-time (continuous) limit only if the damping factor is applied. It now 
depends on a specific quantity how it approaches its infinite-time limiting value. The transition from 
the Fredholm to the Wiener-Hopf solution will then be algebraic for symmetric quantities and weakly 
exponential for antisymmetric functions. Since we are interested rather in the critical exponents than in 
the solution of the Dyson equation (^) itself, we construct long-time limits of the specific combinations 
of the Green functions Ft necessary for the edge behavior. 



1. Orthogonality catastrophe 

To determine the critical edge exponent of the Green function of the core hole we have to evaluate 
the function C(T) from (^). This representation can be rewritten with the aid of the projectors from 
(|4l|). We introduce double indices, one for each independent variable. We can write 

C{T)^- / dA^- ^ e™"'r+'(u;,,u;;;A), (49) 



where wi is either bosonic, Vm, or fermionic, ct;„, frequency according to the subscript Q. Replacing the 
finite-time projectors with their representations using the infinite-time ones (^) we obtain 



C{T) = T / dA- Ve'-"''{r+_(c^„,o.„;A) + e-2''7^F_+(c^„,c^„;A)} 



(50) 



This expression is exact at any time if we insert exact expressions for the functions F-i and F |-. The 

function F |. expresses a contribution from the emission of light, i. e. from filling of the core hole by a 

conduction electron. 

We do not know the exact solution for arbitrary finite times. However, according to the preceding 
section we can use the Wiener-Hopf solution in the long-time regime, i. e. for T ^ T^. from (|4^), to obtain 
an exact long-time asymptotics. In this regime we know the exact asymptotic form of these functions 
to the Wiener-Hopf solution projected onto the finite interval. Although only the symmetric functions 
contribute to the function C{T), we nevertheless have two significant terms in the long-time limit. The 
former one, not damped, is the Wiener-Hopf solution from Section IV. The damped term goes beyond 
the infinite-time Wiener-Hopf solution and makes the difference between the Fredholm and Wiener-Hopf 
solutions. It is at intermediate time scale (T < l/rf) of the same order as the former one and that is why 
the finite-time methods predict a different critical exponent for the Green function of the core hole. 

We can now evaluate nonperturbatively the correction to the Wiener-Hopf critical exponent at in- 
termediate (mesoscopic) times. To this end we have to evaluate the negative function F_ within the 
Wiener-Hopf trick. This function is not a solution of (||) but of a "conjugate" equation where the inte- 
gration interval [0, T] is replaced by [— T, 0]. Then a Wiener-Hopf solution to such an equation is formally 
identical with (|lj), only the exponential functions have inverse meaning. We have to use the negative 
projection 

<i>_(i/,„) = exp j--ln(l + AGe(i.„0) ~ y E , _^ ^ In (l + AG',(u;„) j \ (51) 
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and similarly for the fermionic frequencies. 

We do the same calculation as with the Wiener-Hopf solution in Sec. IV. A. It yields the function 
C(T) in the following form 



CiT) 



iT 



dLL>S{uj) 



27r2 



In^T 



(52) 



This is an exact long-time asymptotic form of the Green function of the core hole at times much greater 
than typical times of the MND Hamiltonian, i. e. ^ 1. Function ( ^2| ) shows a quasicritical behavior 
in the Fredholm regime, i. e. when rjT <gc 1, with a critical exponent of the Nozieres and De Dominicis 
solution. The Wiencr-Hopf critical exponent is realized only for very long times of order of the effective 
lifetime of the transient, excited electron, T ~ r. Nevertheless the Wiener-Hopf solution is able to 
produce an exact quasicritical exponent and nonuniversal prefactor of the core-hole Green function. It 
is necessary to remind that (^2|) was derived exactly only in the weak-coupling regime where no bound 
states interfere. The same applies to the Nozieres and De Dominicis solution as well. 



2. Final-state interaction and absorption amplitude 



A fundamental quantity for determination of the critical behavior of the excited conduction electron 
is the nonequilibrium Green function (^). We can generally write it as 

rT(r,0;ki,k2) = ^5]C^ J2 GcK,ki)r+,K,zi;i;k2). (53) 

We see that this expression does contain antisymmetric combination of the finite-time fermionic and 
bosonic functions. It consists of two contributions. First one is the Wiener-Hopf solution from Sec. IV. 
It contributes due to the difference between arguments of fermionic and bosonic functions, i. e. due to 
the difference w„ — u^. The second one arises from the functional difference between the fermionic and 
bosonic projectors, i. e. due to the difference function AG'r(w) from (^8|). There is no Wiener-Hopf 
solution for this function. We hence cannot derive analogous asymptotic formula for the nonequilibrium 
Green function of the conduction electrons as we did for the core hole in (p^. One has to go beyond 
the infinite-time Wiener-Hopf approach to derive an asymptotically exact behavior of the the function 
Tt(T, 0; ki, k2) in the long-time regime. To this end it is necessary to construct a dynamical equation for 
the difference function Tb — Tf which has the infinite-time limit only with the damping factor e^''-^. Since 
we are unable to find the genuine long-time asymptotics of the Green function of the transient conduction 
electron ( |53| ) at intermediate times, we are also unable to make a conclusion about asymptotic behavior 
of the absorption amplitude for frequencies not too close to the threshold, i. e. for Aw ^ rj. 

We do not have a formula for the absorption amplitude at experimentally relevant frequencies and 
hence we can ask what the physical meaning the Wiener-Hopf solution has. The proper meaning of the 
Wiener-Hopf solution for the absorption amplitude is that it defines the limiting threshold value of the 
edge singularity. It means that it determines the height of the threshold peak at finite temperatures and 
for finite volumes. Behavior of the edge peak observed when limiting the temperature (volume) to zero 
(infinity) is determined by the Wiener-Hopf solution. Also the shift of the ground-state energy is identical 
with the prediction of the Fummi theorem only within the Wiener-Hopf infinite-time solution. Different 
results for different limiting processes and a crossover from intermediate long to very long times give 
explanation to the inconsistency in the interpretation of the critical exponenets from Sec. HI. The limit 
to the infinite volume is at the threshold governed by the Wiener-Hopf solution, while the long-time or 
low-frequency limits are dominated by the quasicritical, finite-time solution of the basic one-body integral 
equation (^. 
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VI. CONCLUSIONS 



In this paper we analyzed in detail the MND problem within the one-body formulation of Nozieres and 
De Dominicis with field-theoretic Green functions and Kubo formula for the absorption amplitude. We 
discussed the differences between two possible ways to solve the fundamental Dyson equation (^) in the 
long-time limit, the finite-time (Fredholm) and infinite-time (Wiener-Hopf) methods. We showed that 
each method is based on a different decomposition of the full function in order to diagonalize (factorize) 
convolution in time. For the finite-time method it is a decomposition into fermionic and bosonic parts, 
the Fourier transforms of which are defined on odd, even frequencies, respectively. A decomposition into 
projections onto retarded (positive-time) and advanced (negative-time) functions is for the Wiener-Hopf 
method appropriate. We determined respective ranges of validity and reliability of both the methods 
and found corresponding small parameters for each of them. We showed that they are in some sense 
complementary. We proved that the solution of Wiener and Hopf for r+ G L J is obtained as a long-time 
limit of the Fredholm solution if ST » 1, where ^ is an effective bandwidth of the conduction electrons. 
It means that the Wiener-Hopf solution is stable to perturbations in the long-time limit and one can 
systematically expand around it. 

We derived the long-time asymptotic limits of all the relevant quantities of the x-ray problem within 
the Wiener-Hopf solution. We calculated the Green function of the core hole, of the transient conduction 
electron excited by absorption of light in the long-time limit and determined not only the critical exponents 
but also the noncritical prefactors. In case of the divergent edge we also obtained the Wiener-Hopf 
expression for the x-ray absorption amplitude. The critical exponents come out in the Wiener-Hopf 
solution as one half of the Nozieres and De Dominicis result. 

To understand the difference we analyzed the role various very large time scales play for the critical 
edge behavior. We showed that it is necessary to introduce an effective lifetime, r, of the transient 
electron-hole pair and compare the relaxation time, T, with this scale in both of the approaches. The 
Fredholm, finite-time approach is characterized by the ratio T/r — > while the Wiener-Hopf solution 
by T/r oo. Unitarity of the theory demands to limit the lifetime r to infinity in the end of the 
calculations. We hence confirmed that the Winer-Hopf results for the relevant quantities of the x-ray 
problem are exact for times comparable with r or frequencies Auj ~ rj = 1/t. We showed that the 
Wiener-Hopf theory is exact for any quantity if the parameter driving the system to the edge singularity 
is of order of the lifetime of the transient electron-hole pair. 

Because of various large scales in the x-ray problem the results we obtain for the critical behavior are 
dependent on the trajectory along which we approach the critical point in the space of relevant parameters. 
Most importantly we found a wide interval of intermediate, experimentally relevant times {T <C r) or 
frequencies [Alu 3> 77) showing a quasicritical behavior. This intermediate-time asymptotics deviates from 
the long-time Wiencr-Hopf one, which causes the differences between the critical exponents calculated 
from the Fredholm or the Wincr-Hopf solutions. We derived general representations for the quasicritical 
long-time behavior of the Green functions of the core hole and of the transient conduction electron. It is 
clear from them that the quasicritical behavior deviates macroscopically from the infinite-time solution. 

Based on the proof of asymptotic exactness of the Wiener-Hopf solution for the functions r± we 
were able to derive an exact long-time asymptotic formula for the Green function of the core hole. This 
formula covers the Wiener-Hopf as well as the quasicritical Fredholm long-time asymptotics. The critical 
exponent for the quasicritical behavior, coincides with the result of Nozieres and De Dominicis in the 
finite-time, Fredholm regime, T/r — > 0. 

With our analysis we demonstrated usefulness and effectivity of the Wiener-Hopf approach to obtain 
the critical edge behavior within the MND model with all its details. We determined the range of 
validity of the Wiener-Hopf solution and explained the differences between its critical exponents and the 
quasicritical exponents from the finite-time approaches at least at weak coupling. Although the Wiener- 
Hopf critical exponents have less experimental relevance, we showed that only the Wiener-Hopf method 
is able to produce an exact, physically relevant long-time asymptotics (exponent and prefactor) of the 
core-hole Green function. We also obtained new results for the nonequilibrium Green function of the 
transient conduction electron and for the absorption amplitude very close to the edge. What remains to 
do is to extend the Wiener-Hopf approach to be able to determine quasicritical behavior of the absorption 
amplitude and to assess effects of bound states to continue the results to strong coupling. 
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